14 research outputs found
Quasi-Monte Carlo rules for numerical integration over the unit sphere
We study numerical integration on the unit sphere using equal weight quadrature rules, where the weights are such
that constant functions are integrated exactly.
The quadrature points are constructed by lifting a -net given in the
unit square to the sphere by means of an area
preserving map. A similar approach has previously been suggested by Cui and
Freeden [SIAM J. Sci. Comput. 18 (1997), no. 2].
We prove three results. The first one is that the construction is (almost)
optimal with respect to discrepancies based on spherical rectangles. Further we
prove that the point set is asymptotically uniformly distributed on
. And finally, we prove an upper bound on the spherical cap
-discrepancy of order (where denotes the
number of points). This slightly improves upon the bound on the spherical cap
-discrepancy of the construction by Lubotzky, Phillips and Sarnak [Comm.
Pure Appl. Math. 39 (1986), 149--186]. Numerical results suggest that the
-nets lifted to the sphere have spherical cap
-discrepancy converging with the optimal order of
Point sets on the sphere with small spherical cap discrepancy
In this paper we study the geometric discrepancy of explicit constructions of
uniformly distributed points on the two-dimensional unit sphere. We show that
the spherical cap discrepancy of random point sets, of spherical digital nets
and of spherical Fibonacci lattices converges with order . Such point
sets are therefore useful for numerical integration and other computational
simulations. The proof uses an area-preserving Lambert map. A detailed analysis
of the level curves and sets of the pre-images of spherical caps under this map
is given